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Abstract: We study the Abelian sandpile model onZ,'^. In d > 3 we prove existence of 
the infinite volume addition operator, almost surely with respect to the infinite volume 
limit n of the uniform measures on recurrent configurations. We prove the existence 
of a Markov process with stationary measure fi, and study ergodic properties of this 
process. The main techniques we use are a connection between the statistics of waves 
and uniform two-component spanning trees and results on the uniform spanning forest 
measure on Z''. 
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1 Introduction 

The Abelian sandpile model (ASM), introduced by Bak, Tang and Wiesenfeld p], has 
been studied extensively in the physics literature, mainly because of its remarkable 
"self-organized" critical state. Many exact results were obtained by Dhar using the 
group structure of the addition operators acting on recurrent configurations intro- 
duced in [4j, see for example [HIE] for reviews. The relation between recurrent config- 
urations and spanning trees, introduced by Majumdar and Dhar [22], has been used 
by Priezzhev to compute the stationary height probabilities of the two-dimensional 
model in the thermodynamic limit [25]. Later on, Ivashkevich, Ktitarev and Priez- 
zhev introduced the concept of "waves" to study the avalanche statistics, and made 
a connection between two-component spanning trees and waves [niHU]. In [2S] this 
connection was used to argue that the critical dimension of the ASM is c? = 4. 

From the mathematical point of view, one is interested in the thermodynamic 
limit, both for the stationary measures and for the dynamics. Recently, in [I] the 
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connection between recurrent configurations and spanning trees, combined with re- 
sults of Pemantle [21j on existence and uniqueness of the uniform spanning forest 
measure on Z'^, has led to the result that the uniform measures /ly on recurrent con- 
figurations in finite volume have a unique thermodynamic (weak) limit fi. (Note: for 
d > 4 the limit was only established for regular volumes such as cubes centered at the 
origin. The extension to arbitrary V is given in the appendix of this paper.) In [TH] 
the existence of a unique limit fi was proved for an infinite tree, and a Markov process 
generated by Poissonian additions to recurrent configurations was constructed. 

A natural continuation of [T] is therefore to investigate the dynamics defined on 
/i-typical configurations. The first question here is to study the addition operators. 
We prove that at least in ci > 3, the addition operators a^, x E Z''- can be defined 
on /x-typical configurations. This turns out to be a rather simple consequence of the 
transience of the simple random walk, and we obtain that the avalanche resulting 
from adding a particle at a given site will be locally finite /i-almost surely (all sites 
topple finitely many times). 

Next, in order to construct a stationary process from the infinite volume addition 
operators, it is crucial that the measure /i is invariant under a^- We show that this is 
the case if the avalanche triggered by adding a particle at x is /i-almost surely finite 
(only finitely many topplings occur). In order to establish almost sure finiteness of 
avalanches, we first prove that the statistics of waves has a bounded density with 
respect to the uniform two-component spanning tree. The final step then is to show 
that the component of the uniform two-component spanning tree corresponding to the 
wave is almost surely finite in the infinite volume limit when d > 3. We deduce this 
from known results on the uniform spanning forest f2M E]. The case d = 2 remains 
an important open question. 

Given existence of a^, and stationarity of fi under its action, we can apply the 
formalism developed in [TH] to construct a stationary process which is informally 
described as follows. Starting from a /i-typical configuration rj, at each site x E Z'^ 
grains are added on the event times of a Poisson process with mean (f{x), where 
ip{x) satisfies the condition 



with G the Green function of simple random walk in Z'^. The condition ensures that 
the number of topplings at caused by additions at all sites has finite expectation at 
any time t > 0. 

In this paper we further study the ergodic properties of the infinite volume process. 
We show that tail triviality of the measure fi implies ergodicity of the process. We 
prove that n has trivial tail in any dimension d > 2. For 2 < d < 4 this is a rather 
straightforward consequence of the fact that the height-configuration is a (non-local) 
coding of the edge configuration of the uniform spanning tree, that is, from the 
spanning tree in infinite volume one can reconstruct the infinite height configuration 
almost surely. This is not the case in d > 4 where we need a separate argument. 
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Our paper is organized as follows. We start with some notation and definitions, 
recalling some basic facts about the ASM. In Sections 3 and 4 we prove existence of 
the addition operator when (i > 3, and show invariance of the measure /x, assuming 
finiteness of avalanches. In Section 5 we prove existence of inverse addition operators. 
Sections 6-8 are devoted to establishing finiteness of avalanches in dimensions ci > 3. 
In Section 6 we make the precise link between avalanches and waves, in Section 7 
we prove that all waves are finite if the uniform two-component spanning tree has 
almost surely a finite component. In Section 8 we prove the required finiteness of the 
component in dimensions c? > 3. Finally, in Sections 9 and 10 we discuss tail triviality 
of the stationary measure, and correspondingly, ergodicity of the stationary process. 

The review papers [121 EH EZ] explain many points that are presented in less detail 
here, and may be useful for the reader. 

2 Notation and definitions 

We consider the Abelian sandpile model, as introduced in [2] and generalized by Dhar 
One starts from a toppling matrix A^.^, indexed by sites in Z'^. In this paper 
A will always be the degree minus the adjacency matrix (in other words, minus the 
discrete lattice Laplacian): 

{2d if X = y, 
-1 if|x-?/| = l, 
otherwise. 

A height configuration is a map rj : Z'^ ^ N = {1,2,...}, and a stable height 
configuration is such that ri{x) < A^^ for all x G Z"^. A site where r]{x) > A^.^. is 
called an unstable site. 

All stable configurations are collected in the set = {1, 2, . . . , 2d}^'' . We endow 
Q with the product topology. For V C Z'^, Qy = {1,2,..., 2d}^ denotes the stable 
configurations in volume V. li rj & Q and W C Z*^, then r]w denotes the restriction of 
T] to the subset W. We also use rjw for the restriction of ?7 G fiy to a subset W <^V. 
Given r] G fly, ^ ^ "we let r^y^yc denote the configuration that agrees with r] 
in V and with ^ in V^. We define the matrix Ay as the finite volume analogon of 
A, indexed now by the sites in V. That is, {Av)xy = A^y, x,y G V. Depending 
on the context, we sometimes interpret Ay as a matrix indexed by Z"^. In that case 
(Ay)^,^ = Aj;yl[x G V]I[y G V], where /[■] denotes an indicator function. 

The toppling of a site x in finite volume V is defined on configurations rj : V —>■ N 

by 

TMiy) = viy) - i^vU (2.1) 

A toppling is called legal if the toppled site was unstable, otherwise it is called illegal. 
The stabilization of an unstable configuration is defined to be the stable result of a 
sequence of legal topplings, i.e., 

Sviv) = Tx„ o o • • • o T^riv), (2-2) 
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where all topplings are legal and such that Sv{ri) is stable. That Sv{ri) is well-defined 
follows from [11123], see also [7j. If is stable, then we define Sv{ri) = rj. The addition 
operators are then defined by 

a^yf] = Sviv + ^x), (2.3) 

where Sx{y) is 1 for y = x and otherwise. As long as we are in finite volume, a^y 
is well-defined and a^^yttyy = ayya^y (Abelian property). 

The dynamics of the finite volume ASM is described as follows: at each discrete 
time step choose at random a site X according to a probability measure p{x) > 
0,x G V, and apply axy to the configuration. After time n, the configuration 
is nr=i ^^i,v ■ ■ ■ (iXn,vV where Xi, . . . , X„ is an i.i.d. sample from p. This gives a 
Markov chain with transition operator 

^/(^) = 5Zp(a;)/Kyr7) (2.4) 

A function f : Q ^ W is called local, if it only depends on finitely many coordi- 
nates, that is, there exists finite V C Z'^, and g : fly ^ such that f{ri) = g{riv), 
rj E Q. Local functions are dense in the space of continuous functions on Q with 
respect to uniform convergence. 

Given a function F{V) defined for all sufficiently large finite volumes in Z*^, and 
taking values in a metric space with metric p, we say that limy F{V) = a, if for 
all e > there exists finite W, such that p{F{V),a) < e whenever V ^ W. For a 
probability measure z/, K^, will denote expectation with respect to u. The boundary of 
V is defined by dV = {y E V : y has a neighbour in V^}, while its exterior boundary 
is defined hj deV = {y E V'^ : y has a neighbour in V}. 



2.1 Recurrent configurations 

A stable configuration rj G is called recurrent, if it is recurrent in the Markov chain 
defined in Section El Equivalently, r] is recurrent, if for any x eV there exists n = nx,n 
such that a^yT] = rj. We denote by TZy the set of recurrent configurations. The 
addition operators a^y restricted to TZy have well-defined inverses a~y, and therefore 
form an Abelian group under composition. From this fact one easily concludes that 
the uniform measure pv on TZy is the unique invariant measure of the Markov chain. 
One can compute the number of recurrent configurations: 

|7^v| = det(Ay), (2.5) 

see Another important identity of [4j is the following. Denote by Nv{x,y,ri) the 
number of legal topplings at y needed to obtain a^rj bom rj + Sx- Then the expectation 
of Nv satisfies 

E^^{Nv{x,y,r])) = Gv{x,y) = {Ay^)xy (Dhar's formula). (2.6) 
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From this and the Markov inequahty, one also obtains Gv{x, y) as an estimate of the 
/iy-probabihty that a site y has to be toppled if one adds at x. We also note that 
for our specific choice of A, Gy is (2(i)~^ times the Green function of simple random 
walk in V killed upon exiting V. 

Recurrent configurations are characterized by the so-called burning algorithm [i]. 
A configuration rj is recurrent if and only if it does not contain a so-called forbidden 
sub-configuration, that is, a subset W <^V such that for all x G W: 

Vix)<~ J2 ^^y (2-7) 

y&W\{x} 

From this explicit characterization, one easily infers a consistency property: if 
r] G TZv and W <i^V then rjw G TZw- This suggests a natural definition of "recurrent 
configurations in infinite volume": we say that r/ G f2 is recurrent, if its restriction to 
any finite V is. We denote this set by TZ: 

n = {r] eQ-.riv eUv ioT all finite V C Z'^}. 

2.2 Infinite volume: basic questions and results 

In studying infinite volume limits of the ASM, the following questions can be ad- 
dressed. In this (non-exhaustive) list, any question can be asked only after a positive 
answer to all previous questions. 

1. Do the measures nv weakly converge to a measure /i? Does fi concentrate on 
the set 717 

2. Is the addition operator defined on /i-a.e. configuration r] E 71, and does it 
leave fi invariant? Does the Abelian property still hold in infinite volume? 

3. Can one define a natural Markov process on 7Z with stationary distribution /i? 

4. Does the stationary Markov process of question 3 have good ergodic properties? 

Question 1 is easily solved for the one-dimensional lattice Z, however, /i is trivial, 
concentrating on the single configuration that is identically 2. Hence no further 
questions on our list are relevant in that case. See [20] for a result on convergence to 
equilibrium in this case. For an infinite regular tree, the first three questions have been 
answered affirmatively and the fourth question remained open [T8]. For dissipative 
models, that is when A^x > 2(i, all four questions are affirmatively answered when 
Aj-j; is sufficiently large 

For Z'^, question 1 is positively answered in any dimension d > 2, using a corre- 
spondence between spanning trees and recurrent configurations and properties of the 
uniform spanning forest on Z*^ pQ. The limiting measure fi is translation invariant. 
The proof of convergence in [T] in the case c? > 4 is restricted to regular volumes, such 
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as a sequence of cubes centered at the origin. In the appendix, we prove convergence 
along an arbitrary sequence of volumes using a random walk result |13j . 

In this paper we study questions 2, 3 and 4 for Z'^, d > 3, and all questions are 
affirmatively answered. 

The main problem is to prove that avalanches are almost surely finite. This is 
done by a decomposition of avalanches into a sequence of waves (cf. [lOl [H]), and 
studying the almost sure finiteness of the waves. The latter can be achieved by a 
two-component spanning tree representation of waves, as introduced in (TUIIII]. We 
then study the uniform two-component spanning tree in infinite volume and prove 
that the component containing the origin is almost surely finite. This turns out to 
be sufficient to ensure finiteness of all waves. 

3 Existence of the addition operator 

In this section we show convergence of the finite volume addition operators to an 
infinite volume addition operator when d > 3. This turns out to be easy, but in 
order to make appropriate use of this infinite volume addition operator, we need to 
establish that /i is invariant under its action, and for the latter we need to show that 
avalanches are finite /i-a.s. 

Let a^y denote the addition operator acting on Qy We define a corresponding 
operator acting on Q using the finite (y)-volume rule, that is, grains falling out of 
V disappear. More precisely, forrjEfl and V 3 x, we define (with slight abuse of 
notation) 

ax,vV = {0'xyVv)Vv'=- (3.1) 

Given rj & Q, call Nv{x, y, rf) the number of topplings caused at y by addition at x in 
7], using the finite (V")-volume rule. Then 

ri + 5.x- ^vNv{x,-,vi) = a^yr], r]eQ,xeV (3.2) 

where Ay is indexed by Z'^. 

We start with the following simple lemma: 

Lemma 3.3. Ny{x,y,r]) is a non- decreasing function of V and depends on t] only 
through rjy . 

Proof. Let V C W. Suppose we add a grain at x in configuration 77. We perform top- 
plings inside V until inside V the configuration is stable, using the finite (ly)-volume 
rule. The result of this procedure is a configuration {axyriv)^^, where possibly ^vnw 
is not stable. In that case we perform all the necessary topplings still needed to sta- 
bilize {axyriv)^v'=nw inside W, using the finite (ly)-volume rule. This can only cause 
potential extra topplings at any site y inside V. □ 

From Lemma 13.31 and by monotone convergence: 

E^{snpNvix,y,7])) = \imE^{Nvix,y,r])). (3.4) 

V y 
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By weak convergence of /iy to /i, and by Dhar's formula (12.61) : 



\imEf,{Nv{x,y,r])) = Urn hm^E^,^ {Nv{x,y,r])) 

< lim Hm E^^^(ArH/(x,i/,r7)) (3.5) 

= \\mGwix,y) = G{x,y), 
w 

where G{x,y) = A~J is {2d)^^ times the Green function of simple random walk in 
Z"^. In the last step we used that d > 3, otherwise Gw{x, y) diverges as t This 
proves that for all x, y G Z'^, N{x,y,ri) = supy Nv{x,y,ri) is /i-a.s. finite. Hence 

/i (Vx, y eZ"^ : N{x, y, r]) < 00) = 1. (3.6) 

Therefore, on the event in (13.61] . we can define 

a^r] = V + Sx- '^N{x,-,v)- (3-7) 

It is easy to see that axf] is stable. This is because axTjiy) is already determined by 
the number of topplings at y and its neighbours, and this is the same as it was in a 
large enough finite volume V . By similar reasons, we also get 

ttxT} = \miaxyr]^ /i-a.s., (3.8) 

where a^y is defined in (13.11) . 

From its definition, one sees that is well behaved with respect to translations. 
Let 6x denote the shift on configuration, that is, OxTjiy) = ri{y — x). Then 

ax = 9xoaoo 9^x, (3.9) 

whenever either side is defined. 

Note that with the above definition of ax, there can be infinite avalanches. How- 
ever, if the volume increases, it cannot happen that the number of topplings at a 
fixed site diverges, and that is the only problem for defining a^. (a problem which 
may arise in d = 2). More precisely, an infinite avalanche that leaves eventually every 
finite region does not pose a problem for defining the addition operator. However, 
as we will see later on, infinite avalanches do cause problems in defining a stationary 
process, at least with our current methods. Hence we need extra arguments to show 
that the total number of topplings is finite /x-a.s. 

We define the avalanche cluster caused by addition at x to be the set 

Cx{v) = {yeZ'':N{x,y,7])>0} (3.10) 

We say that the avalanche at x is finite in 77 if Cx{r)) is a finite set. We say that fi has 
the finite avalanche property, if for all x G Z'^, /i(|Cj;| < 00) = 1. 

In Sections 6-8, we prove the following theorem: 
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Theorem 3.11. Assume d > 3. Then ^ has the finite avalanche property, that is, 
fi{\C^\ < oo) = 1 for all x G Z'^. 

In Section 4, we show that Theorem 13.111 has the following consequence. 

Proposition 3.12. Assume d > 3. Then fi is invariant under the action of a^, 
X G U^, that is, for any fi-integrable function f and for any x G Z*^, 

j f{a^ri)d^i = j fir])dn. (3.13) 

Before moving on to the proofs of Theorem 13.111 and Proposition I3.12[ we prove 
some of their easy consequences. 

Integrating (13.71) over /i and using Proposition 13. 121 we easily obtain the following 
infinite volume analogue of Dhar's formula. 

Proposition 3.14. Assume d > 3. Then 

E^{N{x,y,7])) = G{x,y) (3.15) 

At this point, we cannot compose different a^, since is only defined almost 
surely. 

Proposition 3.16. Assume d > 3. There exists a /i-measure one set f2' C 7^ with 
the following properties. 

(i) For any rj E Vt' and x G IJ^, there exists finite Vx{ri) C TJ^, such that for all 

axV = 0'x,wV- 

(a) For any rj E Q' and x El/ we have a^rj G Q' . 

Consequently, for any rj G Q' , and any xi,...,Xn G Z'^, ax^ax^_^ ■ ■ ■ a^^r] is well- 
defined and all avalanches involved are finite. 

Proof. Define 

VLQ = {r]eTZ: \Cx{r])\ < oo for all x G Z^}. 

By Theorem 13.111 and since fi is concentrated on TZ, we have /i(f2o) = 1. Property (i) 
in the proposition is satisfied for all rj G Qq. For (ii), we need to find a subset of Qq 
invariant under all the a^j's. For n > 1, define inductively the sets 

where here denotes inverse image (not to be confused with the inverse operator 
defined later). Since the Qx are measure preserving, it follows by induction that 
fi{Qn) = 1 for all n. Also, a^. maps f2„ into ^In-i- Therefore, Q' = r\n>o^n satisfies 
both properties stated. □ 
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The following proposition shows that the Abelian property holds in infinite vol- 
ume. 

Proposition 3.17. Assume d > 3. Then 

dxCiyV = CiyaxT], r] G Q'. (3.18) 

Proof. By Proposition 13. 161 for 77 G fi' and for W 3 Vy{ri) U Vx^ayT]) U Vx{ri) U Vy^a^ri) 
we have 

a^ayT] = a^^way^wV = ay,wax,wV = ayaxV- (3.19) 

□ 



4 Invariance of /i under 

In this section, we show that n is invariant under the addition operators, if it has 
the finite avalanche property, that is, we show that Theorem 13.111 implies Proposition 

Proof of Proposition \3.1^ (assuming Theorem \ 3.11\) . It is enough to prove the 
claim for / a local function. In that case, we have 



f{axv)dfi = I fiaxyri)diJ, + ei{V, f) 

f{axyr,)diiw + ei(\/, /) + e^iV, W, f) 
j f{ax,wv)dfiw + e,{V, f) + e2{V, W, f) + e,{V, W, f) 



Here ei and €2 can be made arbitrarily small by (13.81) and by weak convergence. We 
also have 

lesiV, W,f)\< 2\\f\\^fiw{ax,wf ^ a.yf). 
Next, by invariance of fiw under the action of a^^w, 

j fiax,wr])dlJ'W = J fdfxw = J fdfx + e^iW, /). (4.1) 

Here, by weak convergence, €4 can be made arbitrarily small. Therefore, combining 
the estimates, we conclude 



f{axr])dii- j f{ri)dn 



< C limsup limsup /ivK('3'x,Ty/ 7^ o,x,vf)- (4.2) 
V w^v 



Define the avalanche cluster in volume W by 

Cx,w{r]) = {y(^W ■ Nw{x, y, 7]) > 0}, G Vt. 
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Let Df denote the dependence set of the local function /. On the event Cx^wiv)^^^^ = 
we have axyt] = a^^wV- Hence, provided Df ^V, we have 

fiw{ax,wf 7^ < fJ'w{C^,w ndV 0). 

The event on the right hand side is a cylinder event (only depends on heights in V). 
Therefore, the right hand side approaches fi{Cx H dV ^ 0) , as W ^ Z'^ . By Theorem 

Em 

lim/x(C^ ndV ^ 0) = /i(|C^| = oo) = 0, 
which completes the proof. 



5 Inverse addition operators 

In this section we prove that has an inverse defined /i-a.s., provided fi has the finite 
avalanche property. In fact, we show that is the limit of finite volume inverses. 
Let a~Y denote the inverse of Uxy acting on TZy- We define a corresponding operator 
acting on 71, by 

(^xyV = {axyVv)vv-, rjeU. 
This is well-defined, since rjy G IZy- 

Recall that under the finite avalanche property. Proposition 13.161 provided a set 
Vl' of recurrent configurations such that for any rj & Q' and every x E 1^'^, there exists 
a finite set Vx{ri) such that a^r] = axy^(n)V- 

Proposition 5.1. Assume d> 3. There exists a ix-measure one set fl" C fl' with the 
following properties. 

(i) For any Tj G Vt" andx G Z*^ there exists finite V = Vx{ri) such thata~\yrj = a~)^rj 
for allW DV. 

(li) If we define a~^r] = a~^^^^^^r], then a-^a^r] = a^a-'^r] = rj for r] e Q" . 
(Hi) As operators on L2{fi), a* = , that is, the are unitary operators. 
Proof. We first prove that 

hm/i {^WDV: a-^r^ ^ a;;^r/) = 0. (5.2) 

We have 

= A* ^ V : a'l^axT] ^ a-]^axy]^ 

= /i (^3iy D V : a'l^axTj 7^ a^J^a^r/ and "iW D V, ax,wV = ^xyv) ^^'^^ 
+ ey 
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Here we used the invariance of /i under in the first step. The last step follows 
because if a^r] = a^^wV = (^x,vV: then 

a^iv^^xV = (^x!w^x,wV = V = %}/0,-,yr] = %)^a^r]- (5-4) 

As for ey we have 

ev < fi {3W 3 V : a^^wV o.x,vV) (5-5) 
which converges to zero as | Z'^, by the finite avalanche property. 

Since the events in (15.21) are decreasing in V , we get that property (i) of the propo- 
sition is satisfied for /i-a.e. rj. Let us define on this set by setting a~^r] = a~^y ^^^i]. 
Then ( 15.41) shows that has an inverse on a full measure set. By standard argu- 
ments, similar to the one in Proposition I3.16[ we can shrink the set Q' appropriately 
to a set Q" of full measure such that Q" is invariant under and for all x G Z"^. 
Then (i) and (ii) will hold for Q". 

The last statement of the proposition is an obvious consequence of the first two. 

□ 

The above proposition proves that as operators on L2{fi), the generate an 
(Abelian) unitary group, which we denote by G. 

6 Waves and avalanches 

The goal of Sections El [3 and M is to prove Theorem 13.11^ saying that fi has the finite 
avalanche property. 

We will decompose avalanches into so-called waves, that correspond to carrying 
out topplings in a special order. We prove that almost surely, there is a finite number 
of waves, and that all waves are almost surely finite. Without loss of generality, we 
assume that the particle is added at x = 0, the origin. Since the site where we add 
will remain fixed throughout Sections El [7] and [HI henceforth we drop indices referring 
to X from our notation, and simply write C = C{ri) for the avalanche cluster Co(?7). 

We start by recalling the definition of a wave from [ini E]. Consider a finite 
volume W 3 0, and add a grain at site in a stable configuration. If the site becomes 
unstable, then topple it once and topple all other sites that become unstable, except 
0. It is easy to see that in this procedure a site can topple at most once. The set of 
toppled sites is called the first wave. Next, if has to be toppled again, we start a 
second wave, that is we topple once again, and topple all other sites that become 
unstable. We continue as long as needs to be toppled. 

We define awiv) to be the number of waves caused by addition at in the volume 
W. By definition, aw is the number of topplings at 0, that is aw{v) = Nw{0,0,ri). 
For fixed W, let Cw{v) denote the avalanche cluster in volume W. We decompose 
Cw as 

awiv) 

Cw{r]) = U E\^{r]), (6.1) 

i=l 
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where ^^(r]) is the i-th wave in W caused by addition at 0. 

We can define waves in infinite volume as we defined the topphng numbers and 
avalanches in Section [3l by monotonicity in the volume. More precisely, the definition 
is as follows. Fix rj E Q, and assume that is unstable in 77 + 6q. By the argument of 
Lemma [3.31 in Section [3], ^^^(r/) is non-decreasing in W, and therefore we can define 
S^(?7) = Uw'^lviv)- Let rji denote the configuration obtained by toppling every site 
in S^(?7) once, that is r/i = limw Ylxesl^ir])'^^ V- Note that carrying out the first 
wave in the unstable configuration 1] + Sq results in exactly r/i + 5o- By the definition 
of a wave, all sites are stable in rji. If is unstable in rji + 60, we consider 
Since this is again nondecreasing in W, we can define = Uw'^wivi)- This is 

the second wave in infinite volume. Let 772 = limvi/ Y[xe3\^{rii)'^x Vi- Then 172 + 
is the result of carrying out the first two waves in infinite volume on 77 + 6q. Note 
that if |S^(r7)| < 00, we have Ely{rii) = S^(?7) for all large W, and consequently, 
2^(77) = lim.w'^wiv)- We similarly define inductively E^rj) = Uw'^wiVi-i) 
rji = limw Y[x£E]^{ri^.i)'^x: Vi-i- Under the assumption jS-'X?])! < 00, 1 < j < i, we 

also have S*(77) = lim.w'^wii]). For convenience, we define or as the empty 
set, whenever such waves do not exist. 

The easy part in proving finiteness of avalanches is to show that the number of 
waves is finite. Since aw{v) is non-decreasing in W, it has a pointwise limit a{ri), 
and as in (13.51) . 

E^(a) < G(0,0) < 00. (6.2) 

This implies a < 00 /i-a.s. 

In order to prove that C(?7) is finite /i-a.s., we show, by induction on i, that all 
sets S*(?7) are finite /i-a.s. We base the proof on the following proposition, proved in 
Sections [7] and [HI 

Proposition 6.3. Assume d > 3. For i > 1 we have 

limlim sup fiwi^w <^V) = 0. (6.4) 

Proof of Theorem \3 . 1 1\ (assuming Proposition ^. 3^ . We prove by induction on i 
that /i(|S*| < 00) = 1, 2 > 1. To start the induction, note that {S^ C V} is a local 
event, hence weak convergence of to /x and Proposition 16.31 with i = 1 imply that 
/i(|S^| < 00) = L Assume now that /^dS-^l < 00) = 1, 1 < j < i Then 

^V)< 2 ^V',l<]<i)+ /i(S^' 2 V for some l<3 <i). (6.5) 

By the induction hypothesis, the second term on the right hand side can be made 
arbitrarily small by choosing V large. For fixed V' , the event in the first term is a 
local event (only depends on sites in V U deV U U deV). Therefore, the first term 
on the right hand side of (16. 5p equals 

Yim^wi^w 2 ^iv '^V',l<j<i)< limsnp fiwi^w 2 V)- (6-6) 

W w 
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Here the right hand side goes to as f Z'', by Proposition I6.3[ showing that 
/x(|S'| < oo) = 1. 

On the event {a < 00} fl < 00, i > 1}, we can pass to the hmit in fl6.ip and 
obtain the decomposition 

Civ) = U - (^)- (6-7) 

i=l 

It follows that < 00) = 1, which completes the proof of Theorem 13. Ill 

7 Finiteness of waves 

In this section we prove Proposition 16.31 saying that waves are finite. The proof is 
based on a correspondence with two-component spanning trees due to Ivashkevich, 
Ktitarev and Priezzhev [LQ\, which is recalled below. The correspondence allows 
us to use some results on the uniform spanning forest that are stated separately as 
Proposition 17.111 below. The argument is completed with a proof of Proposition 17.111 
in Section [HI 

We now describe the representation of waves as two-component spanning trees 
from P, HO]. Consider a configuration rjw G TZw with r]w{0) = 2(i, and suppose we 
add a particle at 0. Consider the first wave, which is entirely determined by the 
recurrent configuration riw\{o} £ T^wxio}- The result of the first wave on \ {0} is 
given by 

Swiv) = (n^i.H^Vwjwuo} e Tlwxio}- (7.1) 

We associate to any ^ G 7^vk\{o} a tree Tw{^). The tree will represent a wave starting 
at in ^. For the definition of the tree, we use Majumdar and Dhar's tree construction 
[22]. 

Denote by W the graph obtained from Z'^ by identifying all sites in Z"' \ [W \ {0}) 
to a single site 6^^ (removing loops). By [22], there is a one-to-one map between 

recurrent configurations ^ G 7^vi/\{o} and spanning trees of W. The correspondence is 
given by following the spread of an avalanche started at 6^. Initially, each neighbour 
of 6^ receives a number of grains equal to the number of edges connecting it to 5^^, 
which results in every site in Vr\{0} toppling exactly once. The spanning tree records 
the sequence in which topplings have occurred. There is some flexibility in how to 
carry out the topplings (and hence in the correspondence with spanning trees), and 
here we make a specific choice in accordance with [10]. Namely, we first transfer 
grains from 6^ only to the neighbours of 0, and carry out all possible topplings. We 
call this the first phase. When we apply the process to ^ = ?7iy\{o}) the set of sites 
that topple in the first phase is precisely 'E^l{r]) \ {0}. Next we transfer grains from 
6^ to the boundary sites of W, which will cause topplings at all sites that were not 
in the wave; this is the second phase. 

The two phases can alternatively be described via the burning algorithm of Dhar 
[5], which in the above context looks as follows. For convenience, let W denote the 
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graph obtained by identifying all sites in Z'^ \ 14^ to a single site 6w That is, W can 
be obtained from W by identifying and 6w, and calling it 6^. We start with all 
sites of W declared unburnt. At step we burn (the origin). At step t, we 

burn all sites y for which C,{y) > current number of unburnt neighbours of y. (7.2) 

The process stops at some step T = T(^). The sites that burn up to time T is 
precisely the sites toppling in the first phase. We continue by burning 6w in step 
T + 1, and then repeating 07.21) as long as there are unburnt sites. 

Following Majumdar and Dhar's construction [22], we connect with an edge each 
y E W \ {0} burnt at time t to a unique neighbour y' (called the parent of y) burnt at 
time t — 1. This defines a spanning subgraph of W with two tree components, having 
roots and 6w- We denote by Tw{^) the tree component having root 0. Since Tw{^) 
does not contain the vertex 6\y, it can be identified with a subgraph of Z'', and we 
will do so in what follows. (Identifying and 6w merges the two trees into a spanning 
tree of W, yielding the usual spanning tree representation of ^.) 

With slight abuse of language, we refer to the two-component spanning subgraph 
as a two-component spanning tree. By observations made earlier, when ^ = r]w\{o}, 
the vertex set of '7vk(?7h/\{o}) is the first wave Ely{ri). 

We can generalize the above construction to further waves as follows. We define 

Swiv) = (n"j'^\w) ^ '^w\{o}, /c > 1,?7 e TZw (7.3) 

i~o 

If there exists a k-th wave, then its result onW \ {0} is given by (17. 3p . Applying the 
above constructions to ^ = S^^{7]), we obtain that the k-th wave (if there is one) is 
represented by 7^(5*^^(77)). 

If now ^ G TZwxio} is distributed according to ^w\{o}: then Tw{^) is a random 
subtree of Z'^. We will prove that this random tree has a weak limit T, which is 
almost surely finite. But first let us show that this is actually sufficient for finiteness 
of all waves. 

Consider the first wave, and let W ^ V. By construction, '^l^iv) is precisely the 
vertex set of T]y(^riw^^Qy), therefore 

fiwi^wiv) ^V)= iiw{Tw{r,w\{o}) <l V). (7.4) 

Here the right hand side is determined by the distribution of 77vk\{o} under This is 
different from the law of ?7iy\{o} under ixw\{o}^ which is simply the uniform measure on 
'Tlwxio}- It is the latter that we can get information about using the correspondence 
to spanning trees. Indeed, under /ivi/\{o}; the spanning tree corresponding to 77m/\{o} 
is uniformy distributed on the set of spanning trees oiW. In order to translate our 
results back to we show that the former distribution has a bounded density with 
respect to the latter. This will be a consequence of the following lemma. 
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Lemma 7.5. Assume d> 3. There is a constant C{d) > such that 

sup < Cid), (7.6) 

vczd I '<'V I 

where the supremum is over finite sets. 

Proof. By Dhar's formula fl2.5l) . 

I^nwl = det(Av\{o}) = det(A'^) 

where Ay denotes the matrix indexed by sites y E V and defined by {Ay)yz = 
{Av\{o})yz for y,z eV \ {0}, and {A'y)o^ = {Ay)^o = M^)- We have 

Ay + P = A'y 

where P is a matrix for which Py^ = unless y, z & N = {u & 'Z'^ : \u\ < 1}. 
Moreover, maxy^^gv' \Pyz\ ^ 2d — 1. Hence 



|7^y| det(Ay) 

We have {GvP)yz = unless z & N . Therefore 

det(/ + G'vP) = det((/ + GyP)^). (7.7) 

By transience of the simple random walk in ci > 3, we have supy sup^^^ G'y(y, 2;) < 
G(0,0) < cxD, and therefore the determinant of the matrix (/ + GvP)n in (17.71) is 
bounded by a constant depending on d. □ 

We note that an alternative proof of Lemma 17751 can be given based on the following 
idea^ Consider the graph W obtained by adding an extra edge e between and 5w 
in W . Then the ratio in (17.61) can be expressed in terms of the probability that a 
uniformly chosen spanning tree of W contains e. By standard spanning tree results 
[HI Theorem 4.1], the latter is the same as the probability that a random walk in W 
started at first hits 5w through e. 

We continue with the bounded density argument. For any configuration ^ G 
TZw\{o} we have 

/^iy(w\{o} = = |{A; G {1, . . . ,2rf} : {k)Qiw\{o} e 7^vy| . (7.8) 



w\ 



Therefore, 



Mw\m=0 < '^^\w'2rf<C, (7.9) 



/^VF\{0}(W\{0} = \'^w\ 

where, by (17.61) . C > does not depend on ^ or on W . From this estimate, it follows 
that 

^J'w{^w{vw\m) 2 y) 



l^W\{0}iJw{j]W\{Q}) % V) 



< C. (7.10) 
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For a more convenient notation, we let denote the probability measure assign- 
ing equal mass to each spanning tree of W, or alternatively, to each two-component 
spanning tree of W. We can view as a measure on {0, 1}'^'* in a natural way, 
where E*^ is the set of edges of Z*^. By the Majumdar-Dhar correspondence [22], 
corresponds with the measure fiw\{o}, and the law of Tw under fiw\{o} is that of the 
component of under . We keep the notation Tw when referring to . 

We are ready to present the proof of Proposition 16.31 based on the proposition 
below, whose proof is given in Section [HI 

Proposition 7.11. (i) For any d > 1, the limit \imw = u^^^ exists. 

(a) Assume d > 3. Denote the component of under u^^^ by T. Then u^^\\T\ < 
oo) = 1. 



Proof of Proposition \6.3\ (assuming Proposition 7.11 ). By Proposition 17. Ill (i). we 
have 



]im^iw\mirwivw\m) ^v) = ^^uPiTw <iv) = p^'^\t % v). (7.12) 



By Proposition 17.111 (ii), the right hand side of (17. 12p goes to zero as V ] 7/, and 
together with fl7.10l) and fl7.4p . we obtain the i = 1 case of fl6.4p . 

Finiteness of the other waves follows similarly. For k >2 we have by fl7.10p 

<c^^w\m{rw{s^w'v)1V) (7.i3) 

= C iiw\{o}{Jw{rjw\{o}) 2 V"), 

where the last step follows by invariance of iJiw\{o} under Ilj^o '^i=VK\{o}- We have 
already seen in fl7.12p that the right hand side of fl7.13p goes to zero as W,V ^ Z'^, 
which completes the proof of Proposition 16.31 



8 Finiteness of two-component spanning trees 

In this section, we complete the arguments for finiteness of avalanches by proving 
Proposition 17.111 which amounts to showing that the weak limit of Tw as W 1 Z'^ 
is almost surely finite. For this, we briefly review below some results on the uniform 
spanning forest; see [3], |T7] for more background. 

The main statement of Proposition 17.111 is part (ii). It can be deduced from a 
well-known theorem, namely that all trees in the uniform spanning forest in Z*^ have a 
single end. This fact has been known for more general graphs than Z*^, see [3| Theorem 
10.1]. Russell Lyons informed us (private communication) that he, Ben Morris and 
Oded Schramm have independently proved a more general version of statement (ii) 
in the context of giving a new and more widely applicable proof of the single end 
theorem; see [IB] and [ITJ Chapter 9]. 
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For finite ly C Z'^, letj^w denote the probabihty measure assigning equal weight 
to each spanning tree of W. v\y is known as the uniform spanning tree measure in W 
with wired boundary conditions (UST). We use the algorithm below, due to Wilson 
[29] . to analyze random samples from and vw 

Let G* be a finite connected graph. By simple random walk on G we mean the 
random walk which at each step jumps to a random neighbour, chosen uniformly. For 
a path TT = [ttq, . . . , TTm] on G, define the loop-erasure of tt, denoted LE(7r), as the 
path obtained by erasing loops chronologically from vr. 

Wilson's algorithm. Pick a vertex r G G, called the root. Enumerate the vertices 
of G as Let (5'n^)„>i, 1 < i < k he independent simple random walks 

started at xi, . . . , x^, respectively. Let 

T(^) = min{n > : S^^^ = r}, 

and set 

7^^) = LE(SW[0,T(^)]). 
Now recursively define T^^\ i = 2, . . . , k as follows. Let 

T» = min{n > : ^« G U^<j^a^^^, 

and 

7^^) = LE(^(^)[o,r«]). 

(If Xi G Ui<j<i7*^''\ then 7*^*^ is the single point Xj.) Let T = Ui<i<k7^^^- Then T is a 
spanning tree of G and is uniformly distributed |29j . 

Applying the algorithm with G = W and root 6w gives a sample from uw- Sim- 
ilarly, applying the method with G = W and root 6^ we get a sample from z/^''. 

It will be convenient to think of the latter construction also taking place in W, via 
the one-to-one correspondence between the edges of W and W. Note that under this 
correspondence, a path in W that does not use as an internal vertex, maps to a 

path in W. Hence the two-component spanning tree in W can be built from loop- 
erased random walks by regarding {0,5^/} as the "root". In other words, the walks 
attach either to a piece growing from 0, or to a piece growing from 6w, and these two 
growing pieces yield the two components. 

One can extend Wilson's algorithm to infinite graphs G if random walk on G is 
transient [3]. In this case, one chooses the root to be "at infinity", and note that 
loop-erasure makes sense for infinite paths that visit each site finitely many times. 

The measures uw can be reahzed on the same sample space, {0, 1}'^'', as 
introduced earlier. It is well known that has a weak limit u as W ] Z*^, called the 
(wired) uniform spanning forest (USE) on Z'^ [2H E]- When d > 3, the USE can be 
constructed directly by Wilson's method in Z'^, rooted at infinity [HJ Theorem 5.1]. 

We write u for the random set of edges present under z/, that is, we identify 
u G {0, l}'^ with the set of edges e for which uj{e) = 1. This allows us to view u as 
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a (random) subgraph of E'^. We say that an infinite tree T has one end, if there are 
no two disjoint infinite paths in T. It is known that 

z/(all components of uo are infinite trees) = 1, 
z/(each component of uj has one end) = 1; 

see [31 [IE]. 

Proof of Proposition \ 7.11\ Denote the random set of edges present in the two- 
component spanning tree of W by uj^. Let Wn be an increasing sequence of finite 
volumes exhausting Z*^. If i? is a finite set of edges, [31 Corollary 4.3] implies that 
z/|^^(5 C uj\y^) is increasing in n. This is sufficient to imply the weak convergence 

limiy i/^'* = and the limiting spanning forest a; is uniquely determined by the 
conditions 

u^°\BCcu)= limuPjBC cow J, 

n— >oo 

as B varies over finite edge-sets (see the discussion in ^ Section 5]). This proves part 
(i) of the proposition. 

For part (ii), assume d > 3. The configuration under z/^°^ can be constructed by 
Wilson's method directly, by [31 Theorem 5.1]. Since here is part of the boundary, 
the simple random walks in this construction are either killed when they hit the 
component growing from 0, or they attach to a component growing from infinity. 

Assume now that i^'^'^)(|T| = oo) = ci > 0, and we reach a contradiction. We con- 
sider the construction of the configuration under z/^°^ via Wilson's algorithm. Suppose 
that the first random walk, call it S^^\ starts from x 7^ 0. Write x y to denote 
that X and y are in the same component. Then we have 

z/(°)(x ^ 0) = Pr(5(^) does not hit 0) = 1 - ^7^^ ^ 1 as |x| 00. 

G^(U, Uj 

In particular, there exists an x G Z'^, such that 

i^(°)(|r| = 00, X ^ 0) > ci/2. (8.1) 

Fix such an x. Let B{x,n) denote the box of radius n centered at x. Fix hq such 
that E B{x, Uq), and is not a boundary point of B{x, Uq). By inclusion of events, 
( 18. ip implies 

i,(o)(0 ^ 55(x,n), X ^ 0) > ci/2 (8.2) 

for all n > uq. For fixed n > hq, let yi = x, and let y2, ■ ■ ■ ,yK be an enumeration of 
the sites of dB(x, n). We use Wilson algorithm with this enumeration of sites. Let S^^^ 
and T*^*-* denote the i-th. random walk and the corresponding hitting time determined 
by the algorithm. We use these random walks to analyze the configuration under 
both and z/. 

The event on the left hand side of fl8.2p can be recast as 

{T(i) = 00, 3 2 < j < K such that T^-') < 00, s!^^^ = 0}, (8.3) 
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and hence this event has probabihty at least Ci/2. On the above event, there is a 
first index N, 2 < N < K, such that the walk S^^^ hits B^x^rio) at some random 
time (T, where a < T^^\ Let A be the latter event. Since A contains the event in 
(18. 3p . A also has probability at least ci/2. Let p = p{x,no) denote the minimum 
over z G dB{x, uq) of the probability that a random walk started at z hits x before 
without exiting B{x, uq). Clearly, p > 0, and is independent of n. 

Let B denote the subevent of A on which after time a, the walk S^^'^ hits the loop- 
erasure of S^^^ before hitting (and without exiting B{x, no)). We have Pt{B \ A) > p. 
Now we regard the random walks as generating u. By the definition of A^, on the event 
AnB, the hitting times T^^\ . . . , T^^\ have the same values as in the construction for 
u^^\ since the walks do not hit 0. Moreover, on An B, the tree containing x has two 
disjoint paths from dB{x,nQ) to dB{x,n): one is part of the infinite path generated 
by S'^^\ the other part of the path generated by S^^\ Therefore, the probability of 
the existence of two such paths is at least p(ci/2), for all n > no. However, this 
probability should go to zero as n — > oo, because under z/, each tree has one end 
almost surely. This is a contradicion, proving part (ii) of the proposition. 

9 Tail triviality of fi 

In this section we study ergodic properties of /x. In Section [10] we are going to use 
the d > 3 part of the following theorem. 

Theorem 9.1. The measure fi is tail trivial for any d>2. 

Our proof of Theorem 19.11 is divided into two parts. The argument in the case 
2 < (i < 4 is quite simple, and is given in Section [Hill The case > 4 is quite involved, 
and is given in Sections 19.21 and 19. 3[ 

9.1 The case 2 < d < 4 

Proof of Theorem \ 9.1\ [Case 2 < d < 4]. The proof is based on the fact that the 
uniform spanning forest measure u is tail trivial [S, Theorem 8.3]. Let X C {0, l}'^ 
denote the set of spanning trees of Z"^ with one end. Recall the uniform spanning 
forest measure u from Section [HI It was shown by Pemantle [23] that when 2 < d < 4, 
the measure u is concentrated on X. We can regard any uj E X a.s a. tree "rooted at 
infinity", that is, we call x an ancestor of y if and only if x lies on the unique path 
from y to infinity. 

It is shown in pQ that there is a mapping ip : X ^ Q such that n is the image of 
u under ip. Moreover, ip has the following property. Let T^ = T^ioo) denote the tree 
consisting of all ancestors of x and its 2d neighbours in uj. In other words, T^ is the 
union of the paths leading from x and its neighbours to infinity. It follows from the 
results in [1] that rj^ = {ip{uj))x is a function of T^ alone. 

Assume that f{r]) is a bounded tail measurable function. Then for any n, / is 
a function of {rj^ : ||a;||oo > n} only. This means that /(r^) = fijpiuj)) = gioj) is a 
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function of the family {T^i^u) : ||x||oo > n}. Let J^k = o"(co'e : e fl [—k, = 0). For 
1 < k < n consider the event 

En,k= n {T,n[-k,kf = 0}. 

II II oo 

>n 

Observe that En^k ^ ^k, and 5'J[i?„_fc] is jF^-measurable. Using that u has a single 
end u-a.s., it is not hard to check that for any k > 1 

lim u{En,k) = 1- 

n— >oo 

Letting n — oo, this implies that there is an jF^j-measurable function g^, such that 
g = Qk z^-a.s. Since this holds for any A; > 1, tail triviality of u implies that g is 
constant z/-a.s. Letting c denote the constant, this implies 

Kfiv) = c) = ^(/(^H) = c) = 1, 
which completes the proof in the case 2 < d < A. □ 



9.2 Coding of the sandpile in the case d > 4 

The simple proof in Section 19.11 does not work when d > 4, due to the fact that 
the coding of the sandpile configuration by the USF breaks down. Nevertheless, it 
turns out that a coding is possible if we add extra randomness to the USF, namely, 
a random ordering of its components. Due to the presence of this random ordering, 
however, we have not been able to deduce tail triviality of /i directly from tail triviality 
of u, and we need a separate argument. 

We start by recalling results from pLj. Let X denote the set of spanning forests 
of Z'^ with infinitely many components, where each component is infinite and has a 
single end. The USF measure u is concentrated on [3]. Given x G Z'^ and u ^ X, 
let Tx^\uj), . . . ,Tx''\uli) denote the trees consisting of all ancestors of x and its 2d 
neighbours in oo. Here k = kx{uj) > 1. Each Tx^ is a union of infinite paths starting 
at a; or a neighbour of x, and has a unique vertex Vx'^ where these paths "first meet". 
In other words, Vx^ is the first vertex that is common to all of the paths. Let Fx\uj) 
denote the finite tree consisting of all descendants of Vx^ in Tx\uj). Let JF denote the 
collection of finite rooted trees in Z'^. Let denote the set of permutations of the 
symbols {!,..., /}. 

The sandpile height at x is a function of {Fx\uj) , Vx\uj)}i^i and a random cr^. G 
Sfc, in the following sense. 

Lemma 9.2. There are functions ipi : J^^ x S^, / = 1,2,... such that if ax is a 
uniform random element ofT,k, given uj, then 

V. = i^kA{Fi'\ • • • , {Fi'\ vi% (9-3) 
has the distribution of the height variable at x under fi. 
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Proof. This follows from the proofs of Lemma 3 and Theorem 1 in [T] . □ 

Remark 9.4. Here it is convenient to think of ax as a random ordering of those 
components of u that contain at least one neighbour of x. Then one can also view rj^ 
as a function of {Tx'^}\^^ and a^- 

Next we turn to a description of the joint distribution of {r]x}xeAo for finite Aq C 
Z'^. Let A = AqU deAo. Let C'^^\ . . . , C^'^K K = Ka{uj), denote the components 
of the USF intersecting A. Each C^*^ contains a unique vertex v^^ where the paths 
from A n C^*^ to infinity first meet. Let F^'' denote the finite tree consisting in the 
portion of these paths up to v^^ . In other words, F^^ is the union of the paths from 
A n C'^*^ to v^2^ . Each rooted tree {Fx\vx^), x G Aq, 1 < j < kx is a. subtree of 
some F^\ 1 < i < K and the former are determined by the latter. Let a a G '^k be 
uniformly distributed, given u. For each x G Aq, a a induces a permutation in E^^, 
by restriction. Then the lemma below follows from the results in pQ. 

Lemma 9.5. The height configuration in Aq is a function of {{F'f^ ,v'^^2y\f^Y ^'^^ ^k- 
Moreover, the joint distribution of {ax}xeAo is the one induced by oa- 

□ 

From the above, we obtain the following description of /i in terms of the USF 
and a random ordering of its components. Let G be distributed according to 
V. Given lo, we define a random partial ordering on Z"^ in the following way. 
Let C^^\C*^^\ ... be an enumeration of the components of and let \J\, U2, ... be 
i.i.d. random variables, given uj, having the uniform distribution on [0,1]. Define 
the random partial order -<i^ depending on uj and {f/j}j>i by letting x y ii and 
only if X G C*^*-*, y G C^^^ and Ui < Uj. Even though -<t^ is defined for sites, it is 
simply an ordering of the components of uj. The distribution of is in fact uniquely 
characterized by the property that it induces the uniform permutation on any finite 
set of components, and one could define it by this property, without reference to 
the f/'s. This in turn shows that the distribution is independent of the ordering 
(j(2)^ . . . initially chosen. 

Let Q = {0, Ij^'x^" denote the space of binary relations (where for g G Q we 
interpret q{x,y) = 1 as x ^ y, and q{x,y) = otherwise). We denote the joint 
law of (a;, -<) on A" x Q by u. From the couple {uj, -<), we can recover the random 
permutations ax as follows. If v2'\ . . . , vi'''^ are as defined earlier, then 

(axil), ax{k)) = (ji, . . .,jk) if and only if vi''^ ■ ■ ■ -<u. v^J^\ (9.6) 

The discussion above, and Lemma [9.51 easily implies the following lemma. 

Lemma 9.7. Suppose that {uj, -<^) has distribution v . Let r]x be given by ( \9.3\\ . where 
ax is defined by (19. 6p . Then {rjx}x&'^ has distribution ji. In particular, there is a 
v-a.s. defined function : X ^ Q such that fi is the image ofv under ip. 
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□ 



Before we start the argument proper, we need to recall some further terminology 
from [T]. Given finite rooted trees (F, v) = {Fi, Vi)^^ and a finite set of sites A, define 
the events 

D{v) = {vi, . . . ,vk are in distinct components of u}, 
B{F, v) = D{v) n {Ff = Fi, v'i^ =VifoTl<i< K}, 

We also need versions of these events for finite A C Z*^. The wired UST cua in 
volume A can be viewed as the union of one or more components {x,y G A are 
in the same component if they are connected without using the special vertex 6a). 
Let C^^\ . . . ,C^^ be the list of components intersecting A. We define f^^^ and F^\ 

analogously to the infinite volume case, this time using the components C^K Now we 
define 

F>a{v) = {vi, . . . ,vk are in distinct components of ua}, 
Ba{F, v) = Da{v) n = F„ vf^^ = v,hTl<t< K}. 



9.3 Proof in the case d > 4 
9.3.1 Outline of the proof 

Recall that tail triviality is equivalent to the following [H Proposition 7.9]. For any 
cylinder event E' and e > there exists n such that (with Vn = [—n,n]'^ fl Z'^) for 
any event R' G J-'v^ we have 

\fi{E' n R') - i2{E')i2{R') I < e. (9.8) 

Let E = ip~^{E') and R = ip'^iR'), where ip is as in Lemma I^TTl Suppose that E' 
depends on the sites in the finite set Aq, and put A = Aq U deAQ. 

For the proof, we want to show that E and R "decouple" , if n is sufficiently large. 
We try to achieve this by showing that they can be approximated by events that 
depend on portions of u that are spatially separated. The main difficulty is that 
dependence between E and R also exists due to the ordering -<, and it requires work 
to show that the dependence on -< also decouples. Below we give a rough outline of 
strategy for this. 

By Lemma [931 the occurrence or not of the event E is determined by a collection 
of finite tree subgraphs {F^\v^2^) of on, and an ordering of these trees. We get an 
approximation of the event E, if we consider the contribution of only those config- 
urations for which F^^ ^ Vr, 1 < i < K for some large r. Suppose that we have 
also approximated R by an event that depends on the restriction of u; to V^, where 
r < m < n. Condition on the restriction of u to V^, and also on the restriction of -< 
to this portion of u. In particular, for any Wi,W2 G dVm, our conditioning specifies 
whether Wi -< W2 or not. 
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For fixed 1 < i < j < K, let rci and ttj denote the paths in on from v]^ and f ^ , 
respectively, to dVm- Suppose vTj and tcj end in vertices w{i) and w{j), respectively. 
The conditional probability of {f^"* -< f^"*} is determined by the conditional probabil- 
ity of w{i) -< w{j). If r <^ m, then due to fluctuations in the behaviour of the paths 
Hi and Hj, the conditional probability of the events {w{i) -< w{j)} and {w{j) -< w{i)} 
will be approximately equal, and we obtain the desired decoupling. 

In the next section, we specify suitable approximations of E and R. 

9.3.2 Approximating E and R 

We first have a closer look at the event E. We define 

S{F, V, a) = B{F, v) n {v^^i) -< < v^(k)}, 

gE = {iF,v,a):SiF,v,cr)CE}, 
GE{r) = {{F, V, a) eGe: F,CVrioTl<i< K}. 

The event £" is a disjoint union of S{F,v,a) over {F,v,a) G Qe- By Lemma lOTTl we 
have 

^^iE') = ^E)= Yl ^^^iBiF,v)). (9.9) 

We also define an analogue of S" in a finite volume A. Assume that the relation -<g is 
prescribed on the exterior boundary of A. For any realization of the wired UST cua 
there is a unique extension of -<q into A, denoted ^a, where x ^a y if and only if they 
are connected (in ua) to boundary vertices w{x) and w{y) satisfying w{x) w{y). 
Using this extension, we define 

SA{F,v,a) = Ba{F,v) n {v^^i) -<a <a v^(^k)}- 

We let ^A,-<g denote the law of {uj\, -<a) with boundary condition -<g. 
Introduce 

G = G{r) = {F^ C K for 1 < i < K}, 

where we asume that Ao (1 Vr ^ Vn- Now E (1 G is a. disjoint union of the events 
S{F,v,a) over {F,v,a) G ^£;(r). Since A is fixed, there exists tq^e) such that for 
r > ''"o(^) we have i>{G{rY) < e. The event E (1 G{r) will serve as an approximation 
for E. 

Turning to i?, we define 

n = Hn= U U ^^^^^^ 

x&V^ i=l 

The occurrence of R is determined by the collection of edges joining vertices in Ti. 
together with the restriction of -< to 7i. We also introduce for r < m < n and 
C A C the events 

F = F{n, m) = {Hn n Kn = 0} and Fa = F^in, m) = {V^ = A}. 
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Here F is the event that the portion of 00 determining the sandpile configuration in 
does not intersect Vm- The event Rn F will serve as an approximation of R, as 
mentioned in Section 19.3.11 However, we will further decompose F as the disjoint 
union F = IJa-^a- "^^^ reason is that the conditional law of u inside V, given Fa is 
simple: it is z/a (see Section [9.3.31 below). 

The value of m will be chosen large with respect to r. It is easy to see that there 
exists no(m, e), such that if n > no{m,e) then < e. This is because F{n,m) 

is monotone increasing in n, and n^^_,_^F(?T,, m)^ = 0, since each component of the 
USF has a single end. 

For technical reasons, we will in fact need a further subevent of F, on which, given 
the configuration in 7i, with high conditional probability: 

for all x,y E Vr x y implies x ^ y inside Vm- 

Let 

J = {Wx, y E Vj. : if X ^ y then x ^ y inside Vm}- 

There exists mo(r, e), such that if m > mo(r, e), then < esi, where we have set 

£1 = = e/\QE{r)\- Define the event 

Jo = Fn{u{r\LJn„) <ei}, 

where u-h^ denotes the configuration on the set of edges touching Tin- By Markov's 
inequality, 

^{Jq) < i^{F') + v{F n {v{r I un^) >ei])<e + < 2e. 

Summarizing the above, if r > tq, m > mo(r, e) and n > no{m,e), we have 

\fi{E' n R') -u{E nG n Rn Jo))\ < Se- (9.10) 

In the next section, we obtain a decomposition of the event E n G H Rn Jq, that 
allows us to analyze it via Wilson's method. 

9.3.3 Decomposition of E n G D Rn Jq 

We are going to regard the edges of uj being directed towards infinity. By the definition 
of 7i, there are no directed edges from 7i to T>. Therefore, given the restriction of 
uj to 7i, the conditional law of in "D is that of the wired uniform spanning tree in 
V, that is i>£). One can see this by using Wilson's method rooted at infinity to first 
generate Ti and the configuration on Ti, and then the configuration in V. 

Note that the event Fa only depends on the portion of u outside A. We want 
to rewrite the second term on the left hand side of (19.101) by conditioning on Fa, 
the portion of u outside A, and the restriction of ^ to Z'^ \ A. By the previous 
paragraph, the conditional distribution of {u, ^) inside A is given by i^A,^g, where -<q 
is determined by the conditioning. 
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The above implies 



n G n i? n Jo) = ^ / ua,^^{e nG)du. (9.11) 

Since the integration in (19.111) is over a subset of Jq, in what follows, we assume that 
the boundary condition -<q is compatible with the event Jo, in the sense that it arises 
from a configuration belonging to Jo. The expression z/A,-<a(-E fl G) can be further 
decomposed as follows: 

UA,^^iEnG)= Yl ^A,^aiSAiF,v,a)). (9.12) 



In the remainder of the proof our aim is to show that the summand in (19.121) 
is close to i'{B{F,v))/K\, uniformly in A and the boundary condition, if m is large 
enough. In the next section, we formulate precisely the statement we need as Lemma 
19.131 and prove the theorem given Lemma [9. 131 Finally, in Section [9.3.5l we complete 
the argument by proving Lemma [9. 131 



9.3.4 Decoupling lemma and proof of theorem 

Lemma 9.13. There exists a universal constant G and mi{r,e), such that for any 
{F,v,a) G QEij), fn > max{mo(r, e), mi(r, e)}, C A and any boundary condition 
-<Q compatible with Jq, we have 



iyA,^g{SA{F,v,cT)) - 



K\ 



Proof of Theorem \9.1\ [Case d > A] assuming Lemma \9.13i Given £ > let 
f > f^oi^), ^ > max{mo(r, e), mi(r, e)} and n > no{m,e). Then the estimate in 
(Km . formula (El]) and fl9J2D imply 

\l^A,^g{E nG) -u{EnG)\ < Ge. 

Substituting this into (19.111) . and performing the integral and the sum, we get 

\V{E r\GnRr\Jo)-u{En g)V{r n Jo) | < Ge. 

Due to r > tq^s) and n > nQ{m,e), we have u{G'^) < e and ^'(Jo) < 25, which yields 

\^Ji{E' n R') - ii{E')ii{R!)\ < G'e, 
with a universal constant G' . This proves Theorem 19. II in the case d > A. □ 
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9.3.5 Proof of decoupling lemma 



We prove Lemma fQ-lSl by analyzing the event Bx{F, v) in terms of Wilson's algorithm. 
The proof is similar to the proof of Lemma 3 in pLj , however it does not seem possible 
to use that result directly. Before starting the proof proper, we introduce some 
notation. 

Fix (-Fjjfj)^! and a G S^^. Let A = {yi, . . . ,y\A\}- We apply Wilson's method to 
generate part of the wired UST in A with the following enumeration of sites: 

vi, ...,VK,yi,.. .,y\A\- 

Let S^^\ i = 1, ... be independent simple random walks started at fj. Let 7^^*'' 
denote the loop-erasure of S*^*-* up to its first exit time from A. We define a random 
walk event C\ whose occurrence will be equivalent to the occurrence of B\{F,v), by 
Wilson's method. Since the event D/^{v) has to occur, we require that for i = 1, . . . , i^', 
S''-*-' upto its first exit time be disjoint from ^i<j<i'y\^- In addition, the fact that 
B/^{F,v) has to occur, gives conditions on the paths starting at yi, . . . ,y\A\, namely, 

these paths have to realize the events {F^\v^^) = {Fi,Vi), given the paths {7A'*}ili- 
These implicit conditions define C\. More precisely, the loop-erased walk rji started 
at yi has to equal the path in UjFj from yi to {f 1, . . . ,vk}- The loop-erased walk rj2 
started at 2/2 has to equal the path in UjFj from y2 to {f 1, . . . , vk} U rji, and so on. 

We write Pr for probabilities associated with random walk events, and we couple 
the constructions in different volumes by using the same infinite random walks S^^\ 
We also define the random walk event C, corresponding to B{F,v), analogously to 
the finite volume case. 

Proof of Lemma lQ.lSl Let wj^^ denote the first vertex S'*-*-' visits in Z'^ \ A. Then 
we have 

UA,^^^{SAiF,v,a)) = Pr(CA, wi^^'^^ ^ ■■■ ^ H^j'^^^"). (9.15) 

For r < I < m we consider the event Cvi, and write Ci for short. It is not hard to see 
that limA/[CA] = I[C], Pr-a.s., which implies that for / large enough, Pr(C/AC) < ei. 
(Here A denotes symmetric difference.) Hence the difference between the right hand 
side of dinSD and 

Pr(a, Wi^^^^^ ^ ■ ■ ■ ^ Wi^^''^^) (9.16) 

is at most 2ei. Recall that A D Vm, and m > I. By conditioning on the first exit 
points from Vi, f l9.16p can be written as 

Pr(Q) Pr(iy|"(')) ^ ■ • • ^ Wi^^""^^ I Wl'\ . . . , W^""^). (9.17) 

The first factor here differs from Pr(C) = v{B{F ,v)) by at most ei. If m is large 
with respect to Z, the value of the second factor is essentially independent of a. This 
is because the distributions of wf^ and W^l^ given 1^/*^ and W^^'^ (respectively), can 
be made arbitrarily close in total variation distance. This implies that the difference 
between (19.171) and 

Pr(C) Pr(wl') ^ ■ ■ ■ ^ Wf^ I Wl^\ . . . , W^/"^)) 
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is at most ei, if m is large enough, uniformly in A. 

Observe that if wj^;^ ^ W^^^ for some 1 < i < j < K, then the event occurs. 
Since the boundary condition -<q is compatible with Jq, we have 

Ft{Ca, Wji!^ ^ Wi^^ for some I < t < j < K) < UA,^s{r) < ev (9.18) 

It follows that for some universal constant C, if m is large enough 

|Pr(CA, Vrl"^'^^ ^ ■ ■ ■ ^ wjC^""^^) - Pt{C)/K\\ < Cei. 

This proves the lemma. □ 



10 Ergodicity of the stationary process 

Arrived at this point, we can apply the results in [18], and we obtain the following. 
Theorem 10.1. Let ip -.1/ ^ (0, oo) he an addition rate such that 

^ (^(x)G'(0,x) < oo. (10.2) 

Then the following hold. 

1. The closure of the operator on L2{fi) defined on local functions by 

LJ=J2^(^)(^^-^)f (10-3) 

is the generator of a stationary Markov process {rjt : t > 0}. 

2. Let Nf{x) denote Poisson processes with rate (p{x) that are independent (for 
different x). The limit 



Tjt = lim 



r] (10.4) 



n 

exists a.s. with respect to the product of the Poisson process measures on Nf 
with the stationary measure fi on the rj & Q. Moreover, rjt is a cadlag version 
of the process with generator L^. 

Let {rjt : t > 0} be the stationary process with generator = J2x'^(^)(^x ~ 
We recall that a process is called ergodic if every (time-) shift invariant measurable 
set has measure zero or one. For a Markov process, this is equivalent to the following: 
if Stf = f for all t > 0, then / is constant /^-a.s. This in turn is equivalent to the 
statement that Lf = implies / is constant /i-a.s. The tail a-field on Q is defined as 
usual: 

•^oo = n (r{r]{x) : |x| > n} (10.5) 

neN 
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A function / is tail measurable if its value does not change by changing the configu- 
ration in a finite number of sites, that is, if 



fiv) = fi^Wv^) 

for every ^ and V ^ Z'^ finite. 

Theorem 10.6. The stationary process of Theorem \10.1\ is mixing. 

Proof. Recall that G denotes the group generated by the unitary operators on 
L2(/i). Consider the following statements. 

1. The process {rjt : t > 0} is ergodic. 

2. The process {rjt : t > 0} is mixing. 

3. Any G-invariant function is /i-a.s. constant. 

4. is tail trivial. 

Then we have the following implications: 1, 2 and 3 are equivalent and 4 implies 3. 
This will complete the proof, because 4 holds by Theorem 19.11 
It is easy to see that on L2{fi), 

L*=J2^i^)K'-I)- (10.7) 

Hence L and L* commute, that is, L is a normal operator. The equivalence of 1 and 
2 then follows from [28| Lemmas 6 and 7] and an adaptation to continuous time. To 
see the equivalence of 1 and 3: invariance of fi under and implies 

= {L*f\f) 

Hence Lf = is equivalent to / being invariant under all Qx and a~^, and thus 
under the action of G. Finally, to prove the implication 4 ^ 3, we will show that a 
function invariant under the action of G is tail measurable. Suppose / : — M, and 
/ = o-xf = o-x^f /i-a.s. for all x G Z'^. There exists a full measure G-invariant subset 
VLq so that the restriction of / to VLq is G- invariant. If rj and C, are elements of VLq and 
differ in a finite number of coordinates, then 

C=\{ ai.(^)-'?(^)r7 (10.9) 

x&l'^ 

and hence /(?]) = /(C)- This implies that / is /i-a.s. equal to a tail measurable 
function. □ 
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A Appendix 



In this section we show how to extend the argument of [T] in the case d > 4 and prove 
hniA fiA = /i. Using the notation of Section 19.2^ let 

XA,i = distA(fi,5A), i = l,...,K, 
Fa = max iX^i-X^J. 

l<i<j<K ■' 

where distA denotes graph distance in the uniform spanning tree cua- We define the 
random permutation a*^ by the requirement: 

= a if and only if ^a,<t(i) < ■ ■ ■ < X^^aiK), 

where we take a fixed but otherwise arbitrary rule to settle ties. Let K{F) = 
maxi<j<x diam(Fj). The required extension follows once we show the following ana- 
logues of [H Eqns. (18) and (19)]. 

lim/XA {Ba{F,v), Ya < K{F)) = 0, (A.l) 

A 

and 

lim/iA {Ba{F,v), a*j, = a, Ya > K{F)) = -^/i {B{F,v)) . (A.2) 

Most of the argument in [1] does apply to general volumes, and here we detail only 
those points where differences arise. We use the notation introduced in Section 19.3.51 
for Wilson's algorithm. 

We start with the proof of f[0) . Let x,y E Z"^ be fixed, and let and ^(2) 
be independent simple random walks starting at x and y, respectively. Let T^^ and 
T^^ be the first exit times from A for these random walks. The required extension of 



f lA.ip follows from an extension of (27) [IJ , which in turn follows from the statement 

limlimsupPr(l-(5< ^ < 1 + 5) =0. (A.3) 
5^0 A V 7^2) y 



Statement (IA.3P is proved in 

For the extension of (]A.2p . we recall from Section [9.21 the events Ba{F,v) and 



B(F,v) defined for a collection {Fi,Vi)fLi. Let S^''\ i = 1, . . . , K he independent 
random walks started at f j, respectively. Let T^^ be the exit time of 5*^*^ from A. 
Also recall the random walk events Ca and C, and that Cm and Tm are short for Ca 
and t|*'' when A = [— m, m]'^ f] . By the arguments in [T], the required extension of 
( ]A.2I) follows, once we show an extension of (32) [Ij, namely that for any permutation 

lim hmPr [CmX^^^ <■ ■< n^""^) = Pr(C)^. (A.4) 

m^oo A V / A ! 

Observe that Cm and the collection T^^^ = — Tm\ i = 1, . . . , K are conditionally 



independent, given {S^'\Tjf,^)}g^. Therefore, using (1A.3I) . the left hand side of ( 1A.4I) 
equals 

lim lim Pr (C^) Pr fe^i^ < ■ ■ ■ < f^^) . (A.5) 
29 



The second probability approaches 1/K\ for any fixed m, and hence the hmit in (lA.Sp 
equals ViiC) / K\. This completes the proof of the required extension of (]A.2p . 
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